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Babylonian Mathematics, Analog Computation, 
and Lunar Surface Impact Speed 


1 Introduction 


For this project, a combination of ancient mathematics and modern 
electronics will be used to estimate the impact speed of an object 
released from rest at various heights (not to exceed 6 meters) above 
the surface of the moon. 


On page 4, Square roots are displayed in table 1 and impact speeds 
are displayed in table 2. 


2 Mathematical modeling 


First, a brief outline of the Babylonian method to extract the non-negative 
square root of a number using three iterations: 


Let R? = N, such that L < N < H, where L is the perfect square just less than N, 
and H is the perfect square just greater than N. 


Assume that Ri is the square root of L and Ru is the square root of H, then a 
Ru = Ri + 1. 


1st iteration: R1= (Ri+ Ru)/2 = (Ri + Ri + 1)/2 = Ri + 1/2. 
24 iteration: R2 = N/R1=N/(Ri + 1/2). 


34 iteration: R3 = (Ri + R2)/2 = [Ri + 1/2 + N/(Ri + 1/2)]/2. 


From elementary physics, v = V(2gh) (when drag is negligible), where g = 
1.625 m/s? (value at the lunar surface) and h = N (in meters). 


Simplifying, v = V(2 x 1.625 m/s? x N) = 1.80 x Rs. 


3a Computer setup (patch cord version) 


—(Ri +1/2) Ri = (Ri +1/2) 
1/2 


















1/2 


R3 = [(Ru +1/2) + N/(Ru +1/2)]/2 O— (Ri +1/2) + N/(Ri +1/2) 


(Z/T + Tu) /N-= 2a- 





op amp with gain = +1.80 v 


Figure 1: Computer setup 


3b Computer setup (IC/discrete component version) 

















Arrow Color Code: 


Red arrows connect to Vcc. 
Blue arrows connect to VEE. 
Green arrows connect to ground. 





Comments: 


N is from variable DC supply. 





440 k = two 220 k’s in series. 




















Figure 2: Basic breadboard layout 
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4 Results (rounded to two decimal places) 


























N R3 = V(N) R3 = V(N) V(N) 
Babylonian method | Babylonian method | via V key on hand- 
via analog computer | via hand-held held calculator 

calculator using 
formula 

1.50 (RL= 1) 1.26 1.25 1,22 

2.50 (Ri = 1) 1.59 1.58 1.58 

3.50 (Ri = 1) 1.92 1.92 1.87 

4.50 (Ri = 2) 2.15 2.15 2.12 

5.50 (Ri =2) 2.35 2.35 2.35 

6.50 (Ri = 2) 2.56 2.59 2.05 














Table 1: Square root comparisons 





N is replaced with 
h (m) 


Estimated v (m/s) 
at impact 
via analog computer 


v (m/s) at impact 
via physics formula 


|v difference (m/s| 
































1.50 2.26 Peal 0.05 
2.50 2.87 2.85 0.02 
3.50 3.47 2.37 0.10 
4.50 3.89 3.82 0.07 
5.50 4.24 4.23 0.01 
6.50 4.61 4.60 0.01 








Table 2: Impact speed comparisons 
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Figure 3: Basic breadboard layout (top view) 
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Figure 4: Basic breadboard layout (side view) 


